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Abstract 

A FORTRAN77 program is presented that calculates fusion cross sections and mean 
angular momenta of the compound nucleus under the influence of couplings between the 
relative motion and several nuclear collective motions. The no-Coriolis approximation is 
employed to reduce the dimension of coupled- channels equations. The program takes into 
account the effects of non-linear couplings to all orders, which have been shown to play 
an important role in heavy-ion fusion reactions at subbarrier energies. 

PACS numbers: 25.70. Jj, 24.10.Eq 
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PROGRAM SUMMARY 



Title of program: CCFULL 

Catalogue identifier: 17.7 (Experimental Analysis - Fission, Fusion, Heavy- ion) 
Distribution format: ASCII 

Computer for which the program is designed and others on which it has been tested: any 
UNIX work-station or PC. The program has been tested on DEC and DEC-Alpha. 

Operating system or monitor under which the program has been tested: UNIX 

Programming language used: FORTRAN 77 

Keywords: Heavy-ion subbarrier fusion reactions, coupled-channel equations, higher or- 
der coupling, no-Coriolis approximation, incoming wave boundary condition, fusion cross 
section, mean angular momentum, spin distribution, fusion barrier distribution, multi- 
dimensional quantum tunneling 

Nature of physical problem 

It has by now been well established that fusion reactions at energies near and below the 
Coulomb barrier are strongly influenced by couplings of the relative motion of the collid- 
ing nuclei to several nuclear intrinsic motions. Recently, precisely measured fusion cross 
sections have become available for several systems, and a distribution of the Coulomb 
barrier, which is originated from the channel couplings, have been extracted. It has been 
pointed out that the linear coupling approximation, which has often been used in coupled- 
channels calculations, is inadequate in order to analyze such high presicion experimental 
data. The program CCFULL solves the coupled-channels equations to compute fusion 
cross sections and mean angular momenta of compound nucleus, taking into account the 
couplings to all orders. 

Method of solution 

CCFULL directly integrates coupled second order differential equations using the modi- 
fied Numerov method. The incoming wave boundary condition is employed and a barrier 
penetrability is calculated for each partial wave. Nuclear coupling matrix elements are 
evaluated by using the matrix diagonalisation method once the physical space has been 
defined. 

Restrictions on the complexity of the program 

The program is best suited for systems where the number of channels which strongly cou- 
ple to the ground state is relatively small and where multi-nucleon transfer reactions play 
less important role compared with inelastic channels. It also relies on an assumption that 
the fusion process is predominantly governed by quantum tunneling over the Coulomb 
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barrier. This assumption restricts a system which the program can handle to that where 
the sum of the charge of the projectile and the target nuclei Z p + Z^ is larger than around 
12 and the charge product Z v Zt less than around 1800. For most of experimental data 
which were measured to aim to extract fusion barrier distributions, this condition is well 
satisfied. The program also treats a vibrational coupling in the harmonic limit and a 
rotational coupling with a pure rotor. The program can be modified for general couplings 
by directly providing coupling strengths and excitation energies. 

Typical running time 

A few seconds for input provided. The computer time depends strongly upon the number 
of channels to be included. It will considerably increase if one wishes to include a large 
number of channels, as for instance 20. 



3 



LONG WRITE-UP 



1. Introduction 

Fusion is denned as a reaction where two separate nuclei combine together to form a 
composite system. When the incident energy is not so large and the system is not so light, 
the reaction process is predominantly governed by quantum tunneling over the Coulomb 
barrier created by the strong cancellation between the repulsive Coulomb force and the 
attractive nuclear interaction. Extensive experimental as well as theoretical studies have 
revealed that fusion reactions at energies near and below the Coulomb barrier are strongly 
influenced by couplings of the relative motion of the colliding nuclei to several nuclear 
intrinsic motions HI] . Heavy-ion subbarrier fusion reactions thus provide a good opportu- 
nity to address the general problem on quantum tunneling in the presence of couplings, 
which has been a popular subject in the past decade in many branches of physics and 
chemistry. 

Thanks to the recent developments in experimental techniques, fusion cross sections 
can now become measured with high accuracy in small energy intervals. Such high preci- 
sion experimental data have generated a renewed interest in heavy-ion subbarrier fusion 
reactions in recent years 0, [|. For instance, they have enabled a detailed study of the 
effects of couplings on fusion reactions through the so called fusion barrier distribution 
||, U and have thus offered a good opportunity to test any theoretical framework for 
subbarrier fusion reactions. 

Theoretically the standard way to address the effects of the coupling between the 
relative motion and the intrinsic degrees of freedom on fusion is to numerically solve the 
coupled-channels equations, including all the relevant channels. In the past, the coupled- 
channels calculations were often performed using the linear coupling approximation, where 
the coupling potential is expanded in powers of the deformation parameter, keeping only 
the linear term. It has been demonstrated that non-linear couplings significantly affect 
the shape of fusion barrier distributions and thus the linear coupling approximation is 
inadequate in quantitative comparison with the recent high quality data of fusion cross 
sections || 0. The program CCFULL includes the couplings to full order and thus it 
does not introduce the expansion of the coupling potential. Since the dimension of the 
coupled-channels equations with full space is in general too large for practical purposes, 
the program employs the no-Coriolis approximation, which is sometimes referred to as 
the isocentrifugal approximation too, to reduce the dimension || [|. For heavy-ion fusion 



reactions, this approximation has been confirmed to work well The program is 

otherwise exact and take full account of the finite excitation energies of intrinsic motions. 
It includes Coulomb excitations and uses the incoming wave boundary condition inside 
the Coulomb barrier. 

2. Coupled-channels equations 

For heavy-ion fusion reactions, to a good approximation one can replace the angular 
momentum of the relative motion in each channel by the total angular momentum J|| 
El. This approximation, often referred to as no-Coriolis approximation or isocentrifugal 
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approximation, is used in the program. The coupled- channels equations then read 



^ + ^ £ - + ^ , < r » + ^7 i - + e "- £ 



^n{r) + V nm(r)tp m (r) = 0, 



ti 2 d 2 J{J + l)ti 2 ( o), , Z P Z T e 2 

- t c n — hi -y n \i ) -r v nmK' ) VmK' ) ~ u ' 

m. 

1) 

where r is the radial component of the coordinate of the relative motion and /i is the 
reduced mass, respectively. E is the bombarding energy in the center of mass frame and 
e„ is the excitation energy of the n-th channel. V nm are the matrix elements of the coupling 
Hamiltonian, which in the collective model consist of Coulomb and nuclear components. 
These two components are detailed in the following section. is the nuclear potential 
in the entrance channel. In the program, the Woods-Saxon parametrisation 



^> = -IT55(F^ ^„K> + <), ( 2 ) 

is adopted for the nuclear potential 

The coupled-channels equations are solved by imposing the boundary conditions that 
there are only incoming waves at r = r min , and there are only outgoing waves at in- 
finity for all channels except the entrance channel (n—0), which has an incoming wave 
with amplitude one as well. This boundary condition is referred to as the incoming wave 
boundary condition (IWBC) |TjJ, and is valid for heavy-ion reactions, where there is a 
strong absorption inside the Coulomb barrier. The program CCFULL adopts the mini- 
mum position of the Coulomb pocket inside the barrier for r min . Practically the numerical 
solution is matched to a linear combination of incoming and outgoing and Coulomb wave 
functions at finite distance r max beyond which both the nuclear potential and the Coulomb 
coupling are sufficiently small. The boundary conditions are thus expressed as 

ip n (r) -> T n exp(-i f k n (r')dr'j r < r min , (3) 
-> H { -\k n r)5 nt0 + R n H { j +) {k n r) r > r max , (4) 
where 




kn(r) = , £ { E - e n - J(J 2 +^ 2 - V N (r) - - V m (r)\ , (5) 

is the local wave number for the n-th channel and k„ = k n (r = oo). H { j~ ] and H { j +) m eq. 
(|) is the incoming and the outgoing Coulomb functions, respectively. 

In order to ensure that there are only incoming waves at r — > r min , the program 



CCFULL solves the coupled-channels equations outwards from r m , n , first by setting \12 \ 



^n(r min ) = 1, ^ m (r min ) = (m=£n), (6) 

T^n^fflta) = -ik n {r m in), -^ra{r min ) = (m ^ n) . (7) 

Since the first derivative of the wave functions at r min has been explicitly written down 
from eq. (Bl), the wave functions at r = r min + h, h being the radial mesh to integrate 
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the equations, can be determined in the Runge-Kutta method. After the wave functions 
at r = r min + h have been thus obtained, CCFULL solves the coupled-channels equations 



from r 



+ h to r 



in the modified Numerov methods |13]|, since the Runge- 



Kutta method may not be so efficient to solve the second order differential equations. The 
modified Numerov method relates the wave functions at n+i = r min + (i + l)h to those 
at Ti and as 



i+l 




(|l- 4, + ^) 2 - 1 }( 1 -^')* , -( 1 -M^ 1 



where A nm {r) is defined by 



n 



J(J + l)h 2 Z P Z T e 



2/ir 2 



+ 



+ e n -E\5 n 



m V nrn (v) 



(9) 



and ip l are the wave functions at ?v 
Let 

Xnm( r ) be the wave function of the m-ih channel thus obtained, i.e. it is tp m (r) 



which satisfies the boundary conditions eq. 



at r 



At r 



Xr. 



expressed by a superposition of the incoming and outgoing Coulomb waves as 



Xr, 



C nm Hj \k m r) + D nm Hj '(fc m r 



(+)/ 



can be 



(10) 



The coefficients C nm and D nm are determined either by matching the logarithmic deriva- 
tives at r max or by matching the ratio of the wave functions at r max — h to those at r max + h. 
Since the modified Numerov methods does not automatically generate the derivative of 
the wave functions, the latter procedure is more suitable here. The coefficients are then 
obtained as 

(i+l) _ #(+)(*+!) (i-l) 
Jm A.nm Jm A.nm 



and 



Dr 



rr(+)(i-l) rr(-)(*+l) tt(+)(*+1) rr(-)(i-l) 
n Jm n Jm n Jm n Jm 



H (-)(i-l) (i+l) _ rr(-)(i+l) (i-l) 
Jm Aran Jm Aran 



rr(-)(i-l) rr(+)(i+l) 
Jm Jm 



respectively. We have defined = Hj +> (k m ■ (r max + h)), etc. and x % nm = 

Xnm(r max + h), etc. This procedure is repeated for all n and m to determine the ma- 
trices C and D. 

The solution of the coupled-channels equations with the proper boundary conditions 
(§) and (H) is given by a linear combination of Xnm a s 



tr(-)(*+l) rj-(+)(i-l) 
Jm Jm 

r(+)/ 



;n) 



(12) 



n 

This equation satisfies the boundary condition @ at r = r min 



i>m{r ri 



H T nX 



nm\' max J 



^ ] (C nm Hj (k m 7 'max) D nm H j (k m T 



. At r = 7v 



(+)/ 



(13) 

it leads to 
(14) 
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By comparing between eqs. and (fH|), one finds 

The transmission coefficients are then finally obtained by 

r » = ( c_1 L- < 16 > 

For many examples, we are interested only in the inclusive process, where the intrinsic 
degree of freedom emerges in any final state. Taking a summation over all possible intrinsic 
states, the inclusive penetrability is given by 

P j(£) = ^M^l| Tn f. (17) 

The fusion cross section and the mean angular momentum of compound nucleus are then 
calculated by 

a fus (E) = E^) = pE( 2 ^+ 1 )W. (18) 
j o j 



<l> = J2Jaj(E)/J2aj(E), 



^J(2J+l)Pj(E)) /(p£(2J + l)P.7(^ 



(19) 



respectively. In the program CCFULL, the summation over the partial wave is truncated 
at the angular momentum whose contribution to the cross section is less than 10~ 4 times 
total cross section. 

3. Coupling matrix elements 

3.1. Rotational coupling 

In this section, we give explicit expressions for the coupling matrix elements V nm (r) 
in eq. (JH). Let us first consider a rotational coupling in the target nucleus. The nuclear 
coupling Hamiltonian can be generated by changing the target radius in the nuclear 
potential @ to a dynamical operator 

R -> R + 6 = R + (3 2 R T Y 20 + f3 4 R T Y 40 , (20) 

1 in 

where Rt is parametrised as r coup A^ , and /3 2 and /3 4 are the quadrapole and hexadecapole 
deformation parameters of the deformed target nucleus, respectively. The nuclear coupling 
Hamiltonian is thus given by 

V N {r, 6) = ^ R . (21) 

l + exp((r-i? -O)/a) 

We need matrix elements of this coupling Hamiltonian between the \n >= |/0 > and 
\m >= \I'0 > states of the ground rotational band of the target. These can be easily 
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obtained using a matrix algebra ||14|| . In this algebra, one first looks for the eigenvalues 



and eigenvectors of the operator O which satisfies 

6\at >= X a \a > . (22) 

In the program CCFULL, this is done by diagonalising the matrix O, whose elements are 
given by 



O i / 5(2/ + l)(2P+iy // 2 I'V 
° n ' ~ V 4^ mT \ 
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The nuclear coupling matrix elements are then evaluated as 

= <I0\V N (r,6)\r0>-Vf?\r)5 nim , 

= E< /0 I« >< «I J '° > V N(r, A«) - vP(r)5 nim . (24) 

a 

The last term in this equation is included to avoid the double counting of the diagonal 
component. 

For the Coulomb interaction of the deformed target, the program CCFULL includes 
up to the second order with respect to fa and to the first order of fa. Contrary to the 
nuclear couplings, the higher order couplings of the Coulomb interaction have been shown 
to play a rather minor role ||. The matrix elements are then given by 



3Z P Z T R^ 5(21+ 1)(2J'+ 1) 



r 3 V 4-7T 




^ Mijl ^i;;;!, (2B) 

The total coupling matrix element is given by the sum of and V^p. 

3.2. Vibrational coupling 

We next consider a vibrational coupling. Ref. discusses all order nuclear couplings 
for the case where the vibration can be approximated by the harmonic oscillator. In 
realistic case, however, phonon spectra are often truncated at some level, and thus the 
intrinsic motion deviates from the harmonic limit even when the excitation energies are 
equal spaced and/or the electro magnetic transitions do not alter in the linear approxi- 
mation. (See ref. |H| for a discussion on differences between the harmonic oscillator and 
the truncated oscillator, i.e. spin systems.) In such a situation, the matrix formalism 
discussed in the previous section still provides a convenient and powerful technique to 
evaluate the coupling matrix elements 0. For vibrational coupling, the operator O in the 
nuclear coupling Hamiltonian is given by 

d = ^=R T (a[ + a X0 ), (26) 
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where A is the multipolarity of the vibrational mode and a\ (a\o) is the creation (annihila- 
tion) operator of the phonon. The matrix element of this operator between the n-phonon 
state | n > and the m-phonon state \m > is given by 

)• (27) 

V47T 



The rest of the procedure to evaluate the nuclear coupling matrix element is exactly the 
same as the rotational case. The operator O is diagonalised in a physical space and then 
the nuclear coupling matrix elements are calculated according to eq. (p4|) . 

The program CCFULL uses the linear coupling approximation for the Coulomb cou- 
pling of the vibrational degree of freedom. The Coulomb coupling matrix elements are 
thus read 

Vj$(r) = ^^^Z P Z T e^(V^S n , m -i + V^<W0- (28) 



Again the total coupling matrix element is given by the sum of Vjj£j and V^} ■ 

3.3. Transfer coupling 

The program CCFULL includes a pair-transfer coupling between the ground states. 
It uses the macroscopic coupling form factor given by [if 



Ftrans{r)=F t ^- 1 (29) 

where F t is the coupling strength. 

4. Program input and test run 

A description of the format for the input parameters is given in table 1. All parameters 
are entered in free format. The first line contains the parameters specifying the system. 
AP (AT) is the projectile (target) mass and ZP (ZT) is the projectile (target) charge. 
The second line is for the coupling Hamiltonian. RP (RT) is the radius parameter r coup 
of the projectile (target) used in the coupling Hamiltonian. Note that this is in general 
different from the radius parameter used in the nuclear potential (0), which is defined in 
the seventh line. IVIBROTP (IVIBROTT) is an option which specifies the property of 
the intrinsic motion of the projectile (target). If it is set to be —1, the projectile (target) is 
assume to be inert and the fifth (the third and the fourth) line will be ignored. The fusion 
cross sections and the mean angular momentum in the absence of the channel coupling 
can be therefore obtained by setting both the IVIBROTP and the IVIBROTT to be —1. 
When IVIBROTP (IVIBROTT) is set to be zero, the CCFULL assumes that the coupling 
in the projectile (target) is vibrational, while if it is set to be one, the rotational coupling 
is assumed. 

The third line is for detailed information on the target excitation. If IVIBROT is 
zero (i.e., the vibrational coupling), the CCFULL reads OMEGAT, BETAT, LAMBDAT, 
and NPHONONT. OMEGAT is the excitation energy of the single phonon state, BETAT 
is the deformation parameter, and LAMBDAT is the multipolarity of the vibrational 
excitation. NPHONONT is the maximum phonon number to be included. For example, 
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if it is two, up to two phonon states are included in the calculation. The CCFULL 
assumes the harmonic oscillator for a vibrational coupling. The excitation energy of the 
n-phonon state is thus given by n times OMEGAT. Sometimes a user may want to use 
a different value of deformation parameter for the nuclear coupling from that for the 
Coulomb coupling. The CCFULL therefore will ask a user interactively before a run 
whether he/she intends to use a different value of the coupling strength for the nuclear 
coupling. If IVIBROTT is one (i.e., the rotational coupling), the CCFULL reads E2T, 
BETA2T, BETA4T, and NROTT. E2T is the excitation energy of the first 2+ state in the 
ground rotational band of the target nucleus, BETA2T and BETA4T are the quadrapole 
and hexadecapole deformation parameters, respectively. NROTT is the number of levels 
in the rotational band to be included. For instance, if it is 3, the 2 + , 4 + and 6 + states 
are included together with the ground state. The CCFULL assumes a pure rotor for a 
deformed nucleus, and the excitation energy of the J + state is given by 1(1 + l)-E2T/6. 

In many applications, there are two vibrational modes of excitations in the target nu- 
cleus. A typical example is the octupole and quadrapole vibrational excitations in 144 Sm. 
The fourth line is for the second mode of excitation in the target nucleus. The mean- 
ing of OMEGAT2, BETAT2, LAMBDAT2 and NPHONONT2 is the same as OMEGAT, 
BETAT, LAMBDAT and NPHONONT, respectively. The second mode is not included 
when NPHONONT2 is set to be zero. OMEGAT2, BETAT2, and LAMBDAT2 are then 
ignored. When NPHONONT2 is not zero, a user will be asked before a run which of the 
mutual excitation channels he/she intends to include in the calculation. 

The fifth line is the same as the third line, but for the projectile excitations. If there 
exist excitations both in the projectile and the target, the CCFULL takes into account all 
the possible mutual excitation channels between the projectile and the target excitations. 

The sixth line is for the pair transfer coupling. QTRANS is the Q-value for the pair 
transfer channel, while FTR is the coupling strength defined by eq. (|29|). NTRANS is the 
number of the pair transfer channel. In the present version of the CCFULL, NTRANS is 
restricted to be either one or zero. If it is zero, the pair transfer channel is not included 
and QTRANS and FTR are ignored. 

The seventh line is for the nuclear potential in the entrance channel (0). VO is the 
depth parameter of the Woods-Saxon potential, RO is the radius parameter r in eq. (H), 
and AO is the surface diffuseness parameter a. 

EMIN, EMAX, and DE in the next line are the minimum and the maximum value 
of the colliding energy in the center of mass frame and the interval in the energy scale, 
respectively. The CCFULL constructs the distribution of partial cross sections aj as a 
function of J if a single value of the energy is entered, i.e. either when EMIN=EMAX or 
DE=0. 

The accuracy of the calculation is controlled by the matching radius RMAX and 
the mesh for the integration DR in the ninth line. For many application, especially 
for asymmetric system such as 16 + 144 Sm, RMAX=30 fm and DR=0.05 fm provides 
sufficiently accurate results. For heavier systems, such as 64 Ni + 92 Zr, RMAX may have 
to be extended as large as 50 fm. 

The test case shows the fusion cross sections and the mean angular momentum of the 
compound nucleus for the 16 + 144 Sm reaction. The projectile nucleus 16 is assumed 
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to be inert, while the single octupole phonon excitation in 144 Sm is taken into account. 
The transfer channel is not included in this calculation. 
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Table 1: Input to the computer code CCFULL. 



Line 1 


AP, ZP, AT, ZT 


Line 2 


KP, 1V1BKU1P, HI, 1V1BHU1 1 


Line 3 


OMEGAT, BETAT, LAMBDAT, NPHONONT (if IVIBROTT=0) 
E2T, BETA2T, BETA4T, NROTT (if IVIBR0TT=1) 


Line 4 


0MEGAT2, BETAT2, LAMBDAT2, NPH0N0NT2 


Line 5 


OMEGAP, BETAP, LAMBDAP, NPHONONP (if IVIBROTP=0) 
E2P, BETA2P, BETA4P, NROTP (if IVIBR0TP=1) 


Line 6 


NTRANS, QTRANS, FTR 


Line 7 


VO, RO, AO 


Line 8 


EMIN, EMAX, DE 


Line 9 


RMAX, DR 



TEST RUN INPUT 



16. ,8. ,144. ,62. 

1.2,-1,1.06,0 

1.81,0.205,3,1 

1.66,0.11,2,0 

6.13,0.733,3,0 

0,0. ,0.3 

105.1,1.1,0.75 

55. ,72. ,1. 

30,0.05 
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160 + 



TEST RUN OUTPUT 

144Sm Fusion reaction 



Phonon Excitation in the targ. : beta_N= 0.205, beta_C= 0.205, r0= 1.06(fm) 

omega= 1.81(MeV), Lambda= 3, Nph= 1 



Potential parameters: V0= 105.10(MeV), r0= 1.10(fm), a= 0.75(fm) 

Uncoupled barrier: Rb=10.82(fm) , Vb= 61.25(MeV), Curv= 4.25(MeV) 



Ecm (MeV) sigma (mb) <1> 



55, 


.00000 


. 97449E-02 


5 


.87031 


56. 


.00000 


0. 


. 05489 


5. 


. 94333 


57. 


.00000 


0. 


.28583 


6. 


.05134 


58. 


.00000 


1. 


.36500 


6 


.19272 


59. 


.00000 


5. 


.84375 


6. 


.40451 


60. 


.00000 


20. 


.59856 


6 


.86092 


61. 


.00000 


52. 


. 14435 


7. 


.81887 


62. 


.00000 


94. 


. 62477 


9 


.18913 


63. 


,00000 


139. 


. 58988 


10, 


.65032 


64. 


.00000 


185. 


.55960 


11, 


. 98384 


65 


.00000 


234. 


.04527 


13, 


. 13045 


66 


.00000 


283. 


.93527 


14, 


. 18620 


67. 


.00000 


333. 


.26115 


15, 


.21129 


68, 


.00000 


381. 


.21017 


16, 


.20563 


69. 


.00000 


427. 


.61803 


17, 


. 16333 


70. 


.00000 


472. 


.48081 


18, 


.08211 


71. 


.00000 


515, 


. 83672 


18, 


.96273 


72. 


.00000 


557, 


.73621 


19, 


. 80734 
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